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Pentagon and triangle shapes in 28 Si and 12 C are discussed in relation with nuclear density 
wave. In the antisymmetrized molecular dynamics calculations, the IC = 5~ band in 28 Si and 
the K n — 3~ band in 12 C are described by the pentagon and triangle shapes, respectively. These 
negative-parity bands can be interpreted as the parity partners of the K w = + ground bands and 
they are constructed from the parity-asymmetric-intrinsic states. The pentagon and the triangle 
shapes originate in let and 3a cluster structures, respectively. In a mean-field picture, they are 
described also by the static one-dimensional density wave at the edge of the oblate states. In 
analysis with ideal a cluster models using Brink-Bloch cluster wave functions and that with a 
£SJ ' simplified model, we show that the static edge density wave for the pentagon and triangle shapes 

can be understood by spontaneous breaking of axial symmetry, i.e., the instability of the oblate 
states with respect to the edge density wave. The density wave is enhanced in the Z = N nuclei 
' due to the proton-neutron coherent density waves, while it is suppressed in Z =fc N nuclei. 

o 

I. INTRODUCTION 

In light nuclei, some of negative-parity rotational bands with high K quanta are discussed in relation with specific 
symmetry of intrinsic states. One of the famous examples is the 3~ state at 9.64 MeV in 12 C, which has been 
discussed for a long time in connection to an equilateral triangle configuration of 3a cluster structure. The 3 _ state is 
the lowest negative-parity state, and its spin is contradict to the naive expectation from shell model calculations. The 
reason for the low-lying 3~ state is understood by the point group D^h symmetry of the equilateral triangle 3a cluster 
structure [1 3] [see Fig.QJa)], which is characterized by the n-fold symmetry with n = 3 of the intrinsic structure. The 
3~ state is interpreted as the band head of the K 7 ' = 3~ band constructed by the parity and total-angular-momentum 
projection from the D$h symmetry of the intrinsic state. Although the 4~ state in the K 71 = 3~ band has not yet 
confirmed, a possible assignment AT for the level at 13.35 MeV was suggested In microscopic 3a cluster models, 
the K 71 = 3~ band is considered to form a parity doublet with the K v = + ground band [f|. 

In 28 Si, the K 7 ' — 5~ rotational band starting from the 5j~ state at 9.70 MeV was reported in 7-ray measurements by 
Glatz et al. @, and it was discussed with a 7a cluster structure with a pentagon shape 0, I n the 7a cluster model 
with Brink-Bloch (BB) a-cluster wave functions [l|, an oblate solution for negative parity shows the D^h symmetry 
[see Fig. [TJb)]. Since a = 5~ band can be constructed from the D^h symmetry of the intrinsic state, the existing 
K n = 5~ band might be an indirect evidence of the pentagon shape and may be regarded as the parity partner of 
the K 7 ' = + ground state. 

In spite of the reasonable description of the K n = 5~ band in 28 Si with the 7a cluster model, the BB a-cluster 
model is too simple to quantitatively describe low-lying energy spectra of 28 Si @. Moreover, the validity of the ansatz 
that 28 Si consists of seven a clusters is not obvious but it should be checked with frameworks without any cluster 
assumptions because a clusters might be dissociated or melted down due to the spin-orbit force in sd-shell nuclei [![. 

Recently, more sophisticated calculations of 28 Si were performed by one of the authors with antisymmetrized 
molecular dynamics (AMD) [IlllTlJ], which is a framework free from cluster assumptions. The calculations reproduce 
well low-lying positive-parity levels of the oblate ground band and the excited prolate band by incorporating the 
spin-orbit force with a proper strength. Interestingly, the K n = 5~ band is constructed from the oblate state with a 
pentagon shape even though existence of any clusters is not assumed in the calculations fill ]. This means that the 
pentagon shape can be induced by a cluster correlation in the oblate intrinsic state of 28 Si. 

From the viewpoint of the symmetry breaking, the pentagon shape in the intrinsic state is regarded as the sponta- 
neous breaking of the axial symmetry. The surface density is oscillating along the edge of the oblate shape, that is, 
the spontaneous symmetry breaking (SSB) occurs in the rotational invariance around the symmetric axis. In relation 
to the SSB concerning density, this structure is associated with density wave (DW) in a nuclear matter, in which the 
SSB of the translational invariance occurs. The DW in nuclear matter has been discussed for a l ong time (l2rll4l |. and 
it was suggested that the one-dimensional DW could be stable in a low density nuclear matter [l4[. The nonuniform 
nuclear matter with density oscillation has also been investigated with cluster models as the a-cluster matter fl5l - fl7{ . 
In analogy to the nuclear matter DW, the pentagon shape can be interpreted as the static one-dimensional DW at 
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the edge of the oblate state. Our aim is to give description of the pentagon and triangle shapes in terms of the DW 
with the wave number five and three to discuss the relation of the 7 a and 3a cluster structures with the static edge 
DW, i.e., the spontaneous axial symmetry breaking of the oblate states. 

In this paper, we report the AMD calculations of 28 Si while focusing on the pentagon shape in the oblate states. 
In analysis of single-particle wave functions of the obtained AMD wave functions, we show that the pentagon shape 
is expressed by one-particle and one-hole excitations having the wave number five on the oblate state and it can be 
interpreted as the one-dimensional DW. To sec the- instability of the axial symmetry with respect to the pentagon 
shape, analyses of ideal cluster models using BB wave functions are performed. Similarly, focusing on the triangle 
shape of the 3a cluster structure, structures of 12 C are also discussed. We introduce a simplified model for the 
one-dimensional DW in the oblate state by truncating active orbits for particle and hole states, and show that the 
proton-neutron coherent DWs in Z = N nuclei promote the instability of the oblate states with respect to the pentagon 
and triangle shapes. The suppression mechanism of cluster structures in neutron-rich nuclei is discussed from the 
viewpoint of proton DW with excess neutrons. 

The paper is organized as follows: In the next section, we explain the AMD calculations for 28 Si and 12 C. Analyses 
with ideal cluster models using BB a-cluster wave functions are given in Sec. IIII1 and those using extended BB wave 
functions in Sec. IIVI Discussions with a simplified model for the one-dimensional DW are given in Sec. El Finally, in 
Sec. I VII a summary and an outlook are given. 




FIG. 1: The schematic figures for spatial configurations of cluster centers of (a) a triangle structure consisting of three a clusters 
in C and (b) a pentagon structure of seven a clusters in 28 Si. 



II. AMD CALCULATIONS OF 28 SI AND 12 C 



A. Method of AMD calculations 



The AMD method has been applied for various nuclei and succeeded to describe shell-model structures and cluster 
structures of ground and excited states in a light-mass region [l8l - l20l | . Here we briefly describe a simple version of the 
AMD method and its application to 28 Si and 12 C 0, 13 focusing on cluster features of the oblate bands. The details 
of the previous AMD calculations for 28 Si are described in Refs. pUdlj], in which the energy levels of the K n = Of, 
A 71 " = Orj", K n = 3 _ and, A' 71 " = 5~ bands are well reproduced. 

An AMD wave function for an A-nucleon system is given by a Slater determinant of Gaussian wave packets, 
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where the i-th single-particle wave function is written as 

</>Zi(ij) oc exp 



Z^ 



(1) 

(2) 
(3) 



Xi is the spin-isospin function and fixed to be p J, p 4, n 1\ or n \.. The spatial part is represented by complex 
variational parameters, Z X i, Z y i, Z Z j, which indicate the centers of the i-th Gaussian wave packet. The parameter v is 
chosen to be v = 0.15 fm -2 and v = 0.175 fm~ 2 for 28 Si and 12 C so as to minimize the energy of the positive-parity 
state. 

In the AMD model, all the centers {Zi,Z2,--- , Z^} of singlc-nucleon Gaussians are treated independently as 
complex variational parameters. Thus, the AMD method is based completely on single nucleons and therefore it 
is free from such assumptions as cluster existence or axial symmetry Nevertheless, if a cluster structure is favored 
in a system, the cluster structure can be described as an optimum solution of AMD wave functions because BB 



3 



cluster wave functions arc included in the AMD model space. For instance, a cluster formation is expressed by the 
concentration of Gaussian centers for four nucleons, p "f , p J,, n j, and n J, at a certain position. 
By using an effective Hamiltonian, 

i i<j i<j<k 

consisting of kinetic terms and two-body and three-body interaction terms as effective nuclear forces, the energy 
variation is performed within the AMD model space to obtain the optimum solutions, which correspond to the intrinsic 
wave functions for low-lying states. As in Refs. [ToL ITT1 . [Hj], the energy variation is done after parity projection by 
operating (1 ± P r ) on the AMD wave function. After the energy variation for (1 ± P r )$AMD(Z) with respect to Z the 
optimized intrinsic wave functions, $amd(Z' + ') and $amd(Z'~)), are obtained for the positive- and negative-parity 
states, respectively. Then, the total-angular-momentum projection, P^ K , * s operated on the obtained AMD wave 
functions, P^ K (1 ± P r )3>AMD, to calculate expectation values of parity and angular-momentum eigenstates, J^. 

The adopted effective nuclear forces are the same as those in Refs. [lO, [H| with which AMD calculations reproduce 
the energy levels of 28 Si. Namely, the MV1 force (case 1) [2l|, which consists of finite-range two-body and zero-range 
three-body forces, with a parameter set Jb = h = 0, m = 0.62) is used for the central force. As for the spin-orbit force, 
the spin-orbit term of the G3RS force [22[ with the strengths ui = —Un = 2800 McV is adopted. Coulomb force is 
approximated by seven-range Gaussians. 

B. AMD results of 28 Si 

In the present work, we focus only on the oblate rotational bands, = Of , K n = 3j~, and K n = 5j~, though 
the prolate excited = 0^~ band exists in 28 Si [l(J. We adopt the oblate solutions of the intrinsic wave func- 
tions $amd(Z' + -') and $amd(Z'~') obtained by the energy variation after positive- and negative-parity projections. 
| £amd(Z' + ' ) ) and $amd(Z ( - - ' ) ) correspond to the intrinsic states of the lowest positive- and negative-parity bands. 
Density distributions of these AMD wave functions are shown in Fig. [3J Interestingly, the wave function $amd(Z^ + - ) ), 
which corresponds to the = + ground band, shows the pentagon shape due to the 7a-like structure even though a 
clusters are not a priori assumed in the framework. We should comment that a clusters in $amd(Z^ + ^) are somehow 
dissociated due to the spin-orbit force as discussed in Ref. From these two intrinsic states $amd(Z^ + ^) and 

^amdIZ*- - ^), we calculate the states by performing the parity and angular-momentum projections and diagonal- 
izing Hamiltonian and norm matrices with respect to -P//^(l ± P r )$AMD(Z^ + '') and P^xi^ ^ P-) < I ) amd(Z^^). Here, 
states with each parity are described by a linear combination of the parity and angular-momentum eigenstates pro- 
jected from both of ^amd^ -1- - 1 ) and $amd(Z < '~- 1 ). The calculated energy levels of 28 Si are shown in Fig. [3] compared 
with the experimental data of the members in the oblate bands, K n = Of , K v = 3^, and if" = 5^f. The experi- 
mental energy levels are reproduced rather well by the calculations. The calculated in-band E2 transition strengths 
also reproduce well the experimental data (see Table. [I]). The ground band, = 0^, and the lowest negative-parity 
band, = 3^, are dominantly constructed from $amd(Z' + ') and $amd(Z^~- ) ), respectively. The K 7 * = band 
is mainly constructed from $amd(Z^ + ') having the pentagon shape though the = 5^ band member states have 
significant mixing with the K n = 3|f band members. This means that the K v = Of and the K* = bands can 
be interpreted as the parity partners constructed from the parity-asymmetric- intrinsic state $amd(Z' + ') with the 
pentagon shape. 

Next we analyze the single-particle wave functions in the pentagon intrinsic state. The single-particle wave functions 
tpi in Eq. ([2J written by Gaussian wave packets are non-orthogonal to each other. We can make linear transformation 
from the set {fi} to an orthonormal basis set {^} keeping the Slater determinant unchanged except for normalization, 
det{<^} oc det{<^} = $amd(Z). From this orthonormal basis, we construct the Hartree-Fock(HF) single-particles 
{</?f F }, which diagonalize the HF single-particle Hamiltonian as described in Ref. (23|. Analysis of {<£>f F } is helpful 
to discuss intrinsic states in a mean-field picture. 

Among the HF single-particle wave functions in the intrinsic wave function, $amd(Z ( - + ''), we find single-particle 
orbits with pentagon density distributions (see Fig. 0]). The pentagon orbits show the parity asymmetry indicating 
mixing of positive-parity and negative-parity components. We extract each parity component from the pentagon 
orbits, and find that about 5% negative-parity component is mixed in the dominant positive-parity component in 
each orbit. As shown in Fig. 21 both the positive- and negative-parity components show donut shapes in their density 
distributions, and the pentagon orbits can be roughly described by a linear combination c0(o i o,±2) + cV(o,o,=f3) with 
|c'| 2 ~ 0.05 for |c| 2 + |c'| 2 = 1, in terms of harmonic oscillator (H.O.) single-particle orbits labeled by quantum 
numbers (n z , n p , mi) in the cylinder coordinates (see Appendix [X] for the expressions of 0(o,o.±2) an d 0(o,o.±3))- I n 
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FIG. 2: Density distributions of the AMD wave functions for the positive- and negative-parity states in 28 Si and 12 C. 
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FIG. 3: (a) Energy levels calculated with AMD calculations and the experimental levels of the K 71 = 0^, = 3 X , and 
K n = 57 bands in 28 Si. (b) Those of the K n = Of and K n = 3^ bands in 12 C. The experimental data are taken from 
Refs. @, 



the pentagon state of 28 Si, totally eight pentagon orbits are found corresponding to c</>( 0i o,±2) + c '0(o.o,=f3) occupied 
by four species of nucleons, p 'f, p 4-, n t ; and n J,. 



C. AMD results of 12 C 



The calculations of 12 C with the AMD are described in Ref. where the = 0^ and K"* = 'i l bands are 
constructed from triangle states with oblate deformations. In the present work, we use the same effective interactions 
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TABLE I: The calculated and the experimental values of E2 transition strengths in Si. The values in Weisskopf unit, 
W.u -5.05 e 2 fm" 4 are listed. The B(E2) values are calculated with the AMD method. The experimental data are taken from 
Ref. H. 
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as those for 28 Si. 

Both of the intrinsic wave functions $amd(Z^ + ^) and $amd(Z^~- ) ) for positive and negative parity show equilateral 
triangle shapes because of the 3a structure. As seen in density distributions shown in Figs, [^c) and [U^d), the 
development of the 3cv-cluster structure is more remarkable in the negative-parity intrinsic state, $amd(Z^), than 
the positive-parity intrinsic state, $amd(Z' + '). To calculate energy levels, we perform parity and angular-momentum 
projections and obtain the rotational = Of and K n = 37f bands constructed from the triangle intrinsic states, 
$amd(Z' + ') and $amd(Z < - - - 1 ), respectively. If we tolerate the difference of the degree of the cluster development 
between positive- and negative-parity states, the = and K w = 37/ bands are roughly interpreted as the parity 
partners of the parity-asymmetric-intrinsic state with the triangle shape as argued in Ref. [5| . Then, we can say that 
the triangle shape of 12 C has an analogy to the pentagon shape of 28 Si constructing the parity partner K n = 0+ and 
K* = 57/ bands. 

In a similar way to 28 Si, we analyze the HF single-particles {<y9f F } in the 4>amd(Z' + ') and find the parity-mixing 
single-particle orbits with the triangle shape. The density distribution of a ip^ F for the triangle shape is shown in 
Fig-Hl^d), and those of the positive- and negative-parity components extracted from this triangle orbit are shown in 
Figs-Ele) and2Jf). Both the positive- and negative-parity components show donut shape densities, and the triangle 
orbit can be roughly interpreted as a linear combination c0( O ,o,±i) 4- c'</>( ,o,t2) w ith |c'| 2 ~ 0.06 for |c| 2 4- |c'| 2 = 1, in 
terms of H.O. single-particle orbits expressed by cylinder coordinates (see Appendix [5J. In the triangle structure of 
12 C, totally eight triangle orbits are found corresponding to c</>( ,o,±i) 4- c'</>(o,o.=f2) occupied by p 7, p i, n 7, and n \.. 
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FIG. 4: (a) Density of the highest single-particle orbit ip^ F in <I>amd(Z' + ') of 28 Si. (b) and (c) Density for the positive- and 
negative-parity components of the highest orbit. Each component is normalized to be 1 by multiplying n p m with l/n± = 
(ip™\(l ± P r ) 2 \^ F ). (d) Density of the highest single-particle orbit in $amd(Z' + ') for 12 C. (e) and (f) Density for the positive 
and negative-parity components of the highest orbit. Each component is normalized to be 1. 



III. ANALYSIS WITH BRINK-BLOCH a CLUSTER MODELS 



As described in the previous section, the pentagon and triangle shapes are found in the AMD results of 28 Si and 
12 C, respectively. The AMD calculations show that these shapes originate in the 7a- and 3a-cluster features in which 
the a clusters are somehow dissociated because of the spin-orbit force. The fact that the 7a- and 3a-like structures 
are actually formed in the AMD without a priori assuming any clusters indicates that these cluster structures are 
favored in 28 Si and 12 C. 

The pentagon and triangle shapes are interpreted as spontaneous breaking of axial symmetry of the oblate 28 Si and 
12 C. In this section, to understand the mechanism of the symmetry breaking, we investigate properties of ideal 7a- 
and 3a-cluster states by using BB a-cluster wave functions 



A. Brink-Bloch a-cluster wave functions 



BB a-cluster wave functions f° r even-even Z = N nuclei with the mass number A = AX are described by 

A a-cluster wave functions consisting of (Os) 4 a clusters [3, Q. The «-th a cluster is located around a certain position 
Si, and $™ is characterized by a spatial configuration of center positions of Xa clusters, {Si, • • • , Sx}- 

A BB a-cluster wave function for a Xa state can be expressed also by an AMD wave function with a 

specific configuration of Gaussian centers {Z}. When Xi is chosen to be p f, p I, n and n 4 for i = • X}, 
i = {X + 1, • • ■ , 2X}, i = {2X + 1, • ■ ■ , 3A}, and i — {3X + 1, ■ ■ • , 4A}, respectively, and Gaussian centers for four 
nucleons (p t, P I, n n \) are common and real values, = Zi + x = %i+2X = ln+zx = S^/y^ (i = 1 • • ■ X), the 
AMD wave function is equivalent to the corresponding BB a-cluster wave function for X a clusters localizing at the 
positions Si, S2, • • • , S^- 



B. Pentagon 7a BB wave function 



Let us consider the pentagon structure of a 7a system. The a-cluster centers Si of a clusters are taken to have the 
pentagon configuration illustrated in Fig. [TJb) as 




for i = 1, ■ ■ ■ ,5, and 

S, = (0, 0, ±d'Vv) (6) 

for i = 6,7. d is the dimcnsionless pentagon size. Thus, the defined is determined by three parameters is,d,d', 
and hence, we denote the pentagon 7a wave functions by ^^(v,d,d r ). 
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Next we explain the relation between and shell model wave functions by transforming the single-particle wave 
functions of $7^ in the expansion with respect to the pentagon size d. In general, when a-cluster centers {S^} are 
located around the origin, the BB wave function can be connected to a H.O. shell model wave function by using 
invariance of a Slater determinant det{0j(rj)} = no det{<^(rj)} under a linear transformation 4>i(r) — > $(r). Here no 
is a normalization factor. 

For oblate systems, we use H.O. single-particle wave functions 0(n z ,n ,mi) m the expression of cylinder coordinates 
described in Appendix El In the small d! limit, the spatial wave functions for spin-up protons in §^(y, d, d') can be 
transformed to det{</>£(rj)}, which is given by the Taylor expansion with respect to the pentagon size d as follows: 

det{0 Zl ,0z 2 ,--- , 0z 7 } = riodet{<// 1 ,02>--- ,4>'t), 
0'i = 0(o,o,o) + 0(d 2 ), 

02 = 0(1,0,0) +0(d 2 ), 

03 = 0(o,o,+i)+O(d 2 ), 

04 = 0(0,0,-1) + 0(d 2 ), 

05 = 0(0,1,0) + 0(d 2 ), 

06 = 0(0,0,-2) - -^0(o,o,+3) + 0(d 2 ), 

07 = 0(o,o,+2) + -4=0(0,0-3) + 0(d 2 ), (7) 



where n has the order 0(d 9 ,d'). We define 0i i2i3 ,4,5,6,7 = 0(o,o,o), <P (1,0,0), 0(o,o,+i), 0(o,o,-i), i> (0,1,0), 0(o,o,+2), 

0(0,0,-2)- 



In the small d limit, single-particle orbits <\>\ approach <^ and the la wave function becomes equivalent 



to the following 0?iw shell model wave function with the s\s 2 v -p%p%{sd)\{sdfl configuration, 

< B [y, d>,d)-> > ee nt J] det{4>'} 0) X TlT , ^X rrT }, (8) 

TO 

where r = {p,n} and a = {t,l}- In this limit, $7^ has the axial symmetric oblate shape. Here after we consider 
only a small d! limit and investigate properties of (u, d, d') as functions of v and d, $7^ (f, d). In particular, the 
symmetry breaking of the oblate shape caused by the finite d is discussed. 

Le us consider the pentagon shape described by $7 ^ with a finite pentagon size d. As the pentagon size d increases, 
the pentagon shape develops and the axial symmetry breaking of the oblate state enlarges. The leading terms of the 
deviation from are contained in 4>'q and <f>' 7 . The orbits <f>' 6 and (j) 7 are the parity-mixed orbits, and they show 

density with the pentagon shape as expressed in the following explicit form of the density, 

C(r)0 6 (r)=^W0 7 (r) 

1 //A 4 „- r V6 2 f. , dp 



2(tt6 2 ) 3 / 2 



© 4 ^ ,t *( 1 + i^f" s(5 * ) + 0(li2) )- (9) 



The second term cos(5</>) gives the density oscillation with the wave number five along the edge of the oblate shape 
and show the pentagon feature. The orbits (f>g and <f> 7 are given by linear combinations of </>(o,o,=f2) and 0(o.o,±3)- 
Due to the mixing of 0(o,o,±3) m 0(o,o,=f2) °f the amplitude d 2 /6, the density of these orbits changes from the axial 
symmetric density to the oscillating density. This is nothing but the symmetry breaking of the rotational invariance 
around the z axis. If we associate the rotational invariance with the translational invariance of a uniform matter, 
and the z-component of the angular momentum m; with the momentum fc, then we find a good correspondence of 
the 4>q and cf> 7 orbits with the single- particle wave functions of the nuclear matter DW proposed by Overhauser [13]. 
In other words, the pentagon shape can be interpreted as the static DW at the edge of the oblate state. As shown 
below, <f$%(u,d) is expressed by coherent particle-hole configurations from <& 7a . 

We show the particle-hole representation of ^^{v, d) below. We assume that 3>7^f is the Hartree-Fock vacuum 
|0)f, and 0(o,o, ±3)^7-0- and 0(o,o,±2)'^to- are the levels above and below the Fermi level, respectively. We define the 
particle and hole operators as 

"ife.TCT ~ L ±fc,TCT> 

b±q,T<T — C T<J,T-(X, (10) 
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where the labels k = 3 and q = 2 indicate mi for particles and holes. When higher order terms, 0(d 2 ), in the 
single-particle wave functions 0' are ig nored, (v,d) can be approximated to be 

X = to and — \ = r — a . As clearly seen, the product of the particle and hole operators, a± k x b± q _ x , brings quanta 
K = ±5. 



C. Triangle 3a BB wave function 

In a similar way to the pentagon 7a state $7^(1/, d), the equilateral triangle 3a state is related to axial symmetry 
breaking of the oblate state in p-shell, and triangle shape is described by parity-mixed orbits. In the BB a-cluster 
wave function $3^ for the 3a structure, the parameters S, (i = 1, • • • ,3) with a triangle configuration are written as 

'2?r \ . r- . /2tt 



S, = (d^cos(y i) , d^sin(yi) , 0) . (12) 



<&f„ is specified by the parameters v, d as $3^(1/, d). d is the dimensionless triangle size. With a proper transformation 
0z ; (r) — > 0i( r ) an d the Taylor expansion of the transformed single-particle orbits 0J(r) with respect to the triangle 
size d, we can rewrite the spatial wave function for three identical nucleons, 

det{0 Zl ,0z 2 ,0z 3 } = no det^i, 02,03}) ( 13 ) 

with 

01 =01 , 

02 = 02 (O) + ^0(O,O,+2) + 0(d 2 ), 

M> A,'(°) ^ 1 r>M2 

03 = <?3 



2 



0(0,0,-2) + OK), (14) 



where 



»i =0(o,o,o), 
H =0(o,o,-i), 

&3 = 0(0,0,+!)- 



(15) 



In the small d limit, $3^(zy, d) becomes equivalent to the 0?iw shell model wave function, 



The orbits <p f 2 and 0g are the parity-mixed orbits, and their density show a triangle shape with the form 
6 *(r)0 6 (r)=0 7 *(r)0' 7 (r) 

1 /r\ 2 .^L^t^Mnf* 



W>W ^' ^ + ^cos(30 ) + O(^j. (17) 
Similarly to the particle-hole representation of the 7a wave function, $3^(^, d) of the order d can be written in the 



particle-hole representation by using alternative definitions k = 2, q = 1, and |0)f = ^js a (f, d), 

d t , t 



« n ( 1 ! at -*,x & -,,- x ) ( 1 + ) i°> f - ^ 
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D. Development of the 7a and 3a states 

We calculate energies of the 7a pentagon state $7^(1/, d) and the 3a triangle state d) as functions of v and 

d. The energies are evaluated by calculating expectation values of the Hamiltonian H c g given in III Al with respect 
to the intrinsic state ^^(V, oQ, and also the positive-parity projected state (1 + P r )^^(v, d), the negative-parity 
projected state (1 - P r )^^{y,d), and the A'" = 0+ projected one P K=0 (l + P r )^f^(v,d). The parity projection 
corresponds to the restoration of the broken parity symmetry of the intrinsic state, and the K = projection restores 
the broken axial symmetry. The energy minimum state of P K=0 (1 + P r )^^(u, d) may relate to the structure of the 
= 0^ band in 28 Si, while that of (1 - P r )^{v, d) corresponds to the K 7 * = 5^ band because (1 - P r )$f®(v, d) is 
equivalent to K 7 ' = 5~ projected state at least in case of small d. Note that the projected states may contain higher 
correlations beyond a single Slater determinant. 

The contour plots of energy surfaces are shown in Fig. [5] as functions of v and d 2 /6. It is found that the energy 
minimum of the energy surface for ^^(V, d) with no projection locates at d 2 /6 ~ 0.1. The finite pentagon size 
d of the energy minimum indicates the development of pentagon shape, namely, the spontaneous breaking of axial 
symmetry in the intrinsic structure. As seen in the minima of the energy surfaces shown in Fig. [5jd) and (c) , the 
development of the pentagon shape enhances in the K 77 = + projected state, and it is largest in the negative-parity 
projected state. 

The value ci 2 /6 indicates approximately the mixing amplitude of the negative- parity component in the pentagon 
orbits, (p' e and <p' 7 , as given in Eq. ([7]). The value d 2 /6 ~ 0.1 at the energy minimum of the positive-parity projected 
state indicates ~10% mixing which is comparable to 5% mixing of the negative-parity component in the pentagon 
orbits in the AMD wave function $amd(Z + ) for 28 Si. The main reason for the smaller mixing in the AMD result 
than that in the ideal 7a cluster model may be a cluster dissociation effects in the AMD calculations. 

We also calculate energy of the equilateral triangle 3a state $3^(1/, d) as functions of v and d. Energies are 
calculated with respect to the intrinsic state 3>f„ (f, d), and the positive-parity projected state (1 + P r )^f^(v, d), 
the negative-parity projected state (1 - P r )$f^(i/, d), and the K* = 0+ projected one P K=a (l + P r )$f®(v,d). The 
energy minimum state of P K ~°(1 + P r )$>f^(i/, d) describes the structure of the K 71 = 0^ band in 12 C, while that of 
(1 - P r )®f®(v,d) corresponds to the if" = 3^ band. 

The contour plots of the energy surfaces are shown in Fig. [5] as functions of v and e? 2 /4. d 2 /4 is approximately 
the mixing amplitude of </>o,o.=f2 m 0o.o,±i as described in Eq. (|13[) . The finite d 2 /4 value of the energy minimum 
indicates the development of 3a cluster structure. The energy minimum of the energy surface for $>^(i/,d) with no 
projection locates at d 2 /4 ~ 0.1. The cluster development slightly enhances in the K n = + projected state, and it is 
most remarkable in the negative-parity projected state. This is consistent with the AMD calculations of 12 C shown 
in Fig. [2] One of the interesting features of (u, d) is that the energy surface is quite shallow against the large 
triangle size d at v ~ 0.20. This corresponds to three a cluster break up. 



E. Roles of the parity and IC projections in a mean-field picture 

As mentioned before, in the particle-hole representation based on = |0)f, (y,d) can be approximately 

expressed as Eq. (jll[) . At least in the oder d, the negative-parity projected state is nothing but a linear combination 
of lp-lh states, 

(1 - P,)|$?>, d)>«4£ (<>UJ- q ,- x < k ,A«,-x) |0>f- (19) 
^ x 

Thus, (1 — P r )&lja(y, d) is described by the coherent sum of the lp-lh states, a?_ k x ^L g _ x |0)f and a* +k x b\_ q _ x |0)f- It 
indicates that the negative- parity state can be described by the K 71 = 5~ vibration mode on the oblate shape. This is 
an alternative interpretation of the K 71 = 5~ band. However, as already discussed before, since the axial symmetry of 
the oblate state is already broken in the intrinsic state before the negative-parity projection (see Fig. [5|, the if" = 5~ 
band is regarded as the static pentagon "shape" instead of the if" = 5~ vibration on the oblate state. 

The positive-parity projected state corresponds to the mixing of 2p-2h states having K — and K = ±10 into the 
dominant $ 7 ^f"' 1 state. The high K components do not affect to the = + ground band, and actually, they arc 
dropped off in the K = projection. Consequently, the K 7 ' = + projected state contains only the 2p-2h states with 
K = 0, 

p*=°(i+p r )i*£(i/,d)) « io) F - 1 EE flt -^ fot - 9 .-x«U X ' & U-x'i°)F- ( 2 °) 

x x' 
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FIG. 5: (a) Energy expectation values of <^^(y,d) plotted as functions of d 2 /6 and v. (b), (c), and (d) Those of the positive- 



parity state, the negative-parity state, and the K 
d! is taken to be a small value, d! I J~v = 0.1 fm. 



state projected from the intrinsic wave function & 7a {v, d). The parameter 



The contained 2p-2h states are k and — k particle pairs and q and — q hole pairs, and they are associated with Cooper 
pairs in BCS theory {24| . Let us remind the reader that the normal pairings in nuclear systems are considered to 
be neutron-neutron pairing and proton-proton one in the spin 5 = channel. However, the 2p-2h terms in Eq. (|20[) 
have not only spin-zero nn and pp pairs but also spin-one np pairs. Namely, a (S,T) = (1,0) (spin-one isoscalar) 
particle-particle pair and a (S, T) = (1, 0) hole-hole pair couple to be totally S = and T = 0, while a (5, T) = (0, 1) 
(spin-zero isovector) particle-particle pair and a (S, T) = (0, 1) hole-hole pair couple to be S = and T = 0. Because of 
the large number of coherent pairs, the K n — + state projected from the 7a cluster state may gain much correlation 
energy. 

In these analyses, we can say that, in both the negative-parity and K n = + states, the coherent particle-hole 
configurations due to the coherent edge DWs of four kinds \ = p f, p I, n f, and n I play an important role in the 
development of the pentagon shape. We will show the importance of the coherence for the SSB in the later sections. 



IV. EXTENSION OF BB a-CLUSTER MODELS 



As discussed in the previous section, the coherent proton and neutron edge DWs are essential to develop the 
pentagon and triangle shapes. In this section, we investigate the development of the pentagon and triangle shapes 
without the proton-neutron coherence by considering a pentagon or triangle neutron structure with a frozen proton 
structure. For this aim, we extend the BB a-cluster wave functions for the pentagon 7a- and the triangle 3a-cluster 
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FIG. 6: (a) Energy expectation values of $™ (y, d) plotted as functions of d 2 /4 and v. (b), (c), and (d) Those of the positive- 
parity state, the negative-parity state, and the K = state projected from the intrinsic wave function $3,J(^, d). The parameter 
d! is taken to be a small value, d! j ^fv = 0.1 fm. 



states as follows: We assume the pentagon configurations of proton and neutron structures for a Z = N = 14 system 
but take the pentagon size d in Eq. ([5]) independently for protons and neutrons. We take an enough small value of 
the pentagon size d p for protons and vary the size d n for neutrons. Thus, the defined wave function can be written in 
the expansion of the order d n as 

$7 a - n (v, d n ) ~n' Y[ det^f^, ■ ■ • , q>' 7 {0) X pa } 

a 

x Yl det^ ^^, • • • , 4>'!> 0) X na , (j)' & X na , cf>' 7 X na }, 

06 = 0(0,0,-2) - -^=</>(0,0,+3) + 0(d n ), 



i ^™ J. 
(0,0, +2 ) + ^ 



0,0,-3) + 0(4). (21) 



In a similar way, we also assume the triangle configurations of proton and neutron structures for a Z = N = 6 
system by taking the triangle size d in Eq. ((T^|) independently for protons and neutrons. Wc take an enough small 
value of the triangle size d p for protons and vary the size d n for neutrons. Then the wave function can be written in 
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the expansion of the order d n as 

§z a -n{v, d n ) ~ n o ]^det{0'/ O) A'p CT 02 (O) A'p (T 03 (O) A'p (T } 

a 

X JJdet{0fV2^</>3^}, 
a 

03 = 3 (O) -y Vo,-2) +0{dl). (22) 

Moreover, we consider the triangle proton structure in a Z = 6 and TV = 14 system to study the proton edge DW 
in a neutron-rich system. For the frozen neutron structure, we adopt a pentagon configuration of the neutron part 
with an enough small pentagon size d n . The proton structure is assumed to be a triangle structure with the triangle 
size dp, which is a variational parameter. The wave function can be written in the expansion of the order d p as 

$2o C _ p (V, dp) w n' det{0' 1 (o) X pa , 4>' 2 X P<T , 4>' 3 X pa } 

x J] det{<j>'t 0) X na ,--- ,0 7 (O) *„ CT }, 

<X={t4} 

02=02 (O) + ^0(O,O, + 2)+O(^), 

03 =03 (O) - f 0(0,0,-2) +0(4). (23) 

This model corresponds to a 3a core structure in 20 C system. 

We calculate energies of $7 Q _ n (^, d n ), <&3 a - n (^, and <&2oc_ p (i/, d p ) states and compare the results with d) 
and (i/, d) . The energies are evaluated by calculating expectation values of the effective Hamiltonian H e s for these 
states with no projection, the positive- and negative-parity projected states, and the K n = + projected states. 

At first, we compare the pentagon size dependence of the energies of <f>7 a . n (v, d n ) having the frozen proton structure 
with that of $y^(^, d) having the proton-neutron coherent pentagon shapes. The energy curves are shown in Fig. [7] 
In each system, v is fixed to be the optimum value at the energy minimum solution in the v-d plane for the positive- 
parity projected state. As already discussed in the previous section, ^^{v, d) shows the deep energy pocket around 
the energy minimum at the finite d value [see Fig. [7ta)] . This indicates the development of the pentagon shape, which 

corresponds to the spontaneous breaking of the axial symmetry of the oblate state The potential pockets 

are deeper in the projected states than the intrinsic state with no projection. In contrast to the energy curve for 
d), the energy curve for $7 Q _„(^, d n ) with no projection has the minimum around d n = 0, which corresponds 
to the axial symmetric oblate state • Even in the projected states, there is no deep pocket in a finite d n region, 

and the pentagon shape of the neutron structure is suppressed in the frozen proton structure. This means that the 
neutron edge DW on the oblate state 3> 7 '^"' > docs not occur without the coherent proton edge DW. It may lead to 
suppression of pentagon shape in Z ^ N nuclei. 

Next we discuss the triangle structures in Z = 6 nuclei. In a similar way to the pentagon structure, we compare 
triangle size dependence of the energies of $3 Q! -n(i / , d n ) having the frozen proton structure with that of $3^ (z/, d) 
having the proton- neutron coherent triangle shapes in Figs. EJa) and[5Jb). Again we find that each energy curve 
for $3 Q -n(^,rfn) has the minimum around d n = 0, which corresponds to . It is in contrast to the features of 

$3^(1/, d), which shows the deep energy pocket at the finite d indicating to the developed triangle shape. We also 
consider the proton triangle shape in the Z = 6, N = 14 system, $2o c _ p (^, d p ). The finite d p corresponds to the 
development of the 3a core structure in 20 C system with the oblate proton and neutron structures. The energy of 
$2oc_p(^, d p ) is the smallest around d = and increases as the triangle size d p becomes large. Compared with the 
energy curve of $3 Q -ri(j / , d n ), the triangle proton structure is significantly unfavored in the Z = 6, N — 14 system. 

Then we conclude that the proton-neutron coherence is essential in development of the pentagon and triangle 
structures of the oblate states in Z = N nuclei. Needless to say, this is consistent with the cluster aspect of Z = N 
nuclei. In the oblate state of neutron-rich C, the triangle cluster structure is suppressed. The first reason for the 
quenching of cluster structure is lack of the proton- neutron coherence. The second reason is the expanded level spacing 
of proton orbits in neutron-rich nuclei because protons are deeply bound due to excess neutrons. Since the energy cost 
for lp-lh proton excitations increases in the neutron-rich system, the correlation energy due to the triangle structure 
may not be able to overcome the cost. We shall discuss the details in the next section. 
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FIG. 7: (a) The energy of the 7a state d) as a function of d 2 /6. d is the pentagon size for protons and neutrons, (b) 

The energy of the Z = N = 14 state <§>7 a -n(y, d n ) with the frozen proton structure as a function of <f 2 /6. d n is the pentagon 
size for neutrons. In each system, v is fixed to be the optimum value at the energy minimum solution in the v-d plane for the 
positive-parity projected states as (a) v = 0.15 fm -2 and (b) v = 0.135 fm -2 . The parameter d' is chosen to be d! /y/v = 0.1 
fm. 



V. EDGE DENSITY WAVE AND SPONTANEOUS SYMMETRY BREAKING 

As already mentioned, the pentagon and triangle structures can be interpreted as the static edge DWs at the 
surface of the oblate states, which is expected to connect with the spontaneous symmetry breaking (SSB) of rotational 
invariance around the symmetric axis. 

In systems with strong interaction, we know various SSB phenomena such as nuclear BCS, chiral symmetry breaking, 
and color superconductivities. These SSB occur both homogeneously and inhomogcncously. In particular, inhomoge- 
neous SSB phases are discussed in the framework of nuclear DWs, chiral DWs, and Fuldc-Fcrrcl-Larkin-Ovchinnikov 
state in color superconducting phase (l2l . [2oT - l39| , whose phase breaks translational invariance and the corresponding 
condensation operator depends on the spatial coordinates. In more general, the SSB resulting in a spatially nonuni- 
form vacuum relates to condensation operators with finite momenta. When we understand nuclear matter DW as the 
instability of the Fermi surface, the condensation operator is given in the form a\ F b kF , which has the momentum 2fcp 
(k,F is the Fermi momentum). In condensed matter physics, inhomogeneous phases with DW arc discussed as charge 
DWs and spin DWs jSEl. 

These phenomena in infinite systems are discussed in terms of the order parameters, which are characterized by 
non-zero expectation values of certain operators. In finite systems, however, the symmetry cannot be broken in the 
energy cigenstates, because the symmetry is restored even if it is broken in the intrinsic state. Nevertheless, it is useful 
to discuss the SSB in the state before projection or restoration by analyzing expectation values of specific operators 
resemble to the condensation operators as done for the BCS phenomena in finite nuclei. For the pentagon and triangle 
structures, the expressions in Eqs. (|TT|) and (fT5|) arc the form similar to the matter DW operators a k fct • 

In this section, we describe the SSB for edge DWs by introducing a simplified model in Appendix [C] and discuss 
the development and suppression of the pentagon and triangle structures from the viewpoint of the edge DW. In 
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FIG. 8: (a) The energy of the 3a state d) as a function of d 2 /4, where d is the pentagon size for protons and neutrons, 

(b) The energy of the Z — N — 6 state <&3 a -n(y, d„) with the frozen proton structure as a function of d 2 /4. d n is the pentagon 
size for neutrons, (c) The energy of the Z — 6, N = 14 system <3>2o c _ p (i/, d p ) with the frozen proton structure. The pentagon 
size d n for the frozen neutron structure is taken to be d 2 — 0.025. In each system, v is fixed to be the optimum value at the 
energy minimum solution in the v-d plane for the positive-parity projected states as (a) v — 0.175 fm~ 2 , (b) v = 0.160 fm~ 2 , 
and (c) v = 0.135 fm~ 2 . 



this model, is assumed to be the Hartree-Fock vacuum |0)f, and the orbits \(f>(o,o,±q)X Ta ) and \(f>(o,o,±k)X T a) 

are considered to be active Hartree-Fock single-particle states. It means that the model space is truncated within 
4>(p,o,±k) %to f° r particle states and 4> '(ofi,±q)X T <7 for hole states. This is equivalent to the model of 8 particles for 16 
states, which is a kind of half filled models. As for the residual interaction, we assume a contact interaction and adopt 
-Hdw defined in Eq. (|C7[) . Note that this model is applicable also to the 3a oblate state by replacing q = 2 and k = 3 
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for the 7a state with q = 1 and k — 2. Then, the Hamiltonian in the particle- hole representation can be written as 

H = Hq + Hi + -ff dw ! 
H = F <0|#|0) F , 



H 1 = y~" E ktT a^ +k x a +k , x + ^ E k , T a*_ 



k, x 



a-h 



E Q,rbi q , x b +q , x ~ 51 E k,rbl q , x b- q , : 



X X 
^DW - 2 2^ G X,X' 

x,x' 

+ a -fc, x 6 -g,- X fe -9-x'a-fc,x'] ■ (24) 

Here x = TCT an d — X = t — c. 

We use an ansatz for the new vacuum of the edge DWs with the axial-symmetry breaking as 

l*> =U{ v - + u Akji q ,- x )U{<-< a Ux b U-x)\°)^ ( 25 ) 
x x 

with 

\v T \ 2 + \u T \ 2 = 1, (26) 

where u r and u T are variational parameters determined by the energy variation, and time reversal invariance is taken 
into account. Our ansatz Eq. (|25j) has the same form as that obtained in the approximation that the quantum 
fluctuation of the particle-hole operators such as a^ fe x b± q _ x are omitted as shown in Appendix |D| It is clear that 

SSf is equivalent to d) with u T = —d/\Q in the order d approximation given by Eq. dHJ) (or <f>f*(v,d) with 

u T = d/2). In general, the coefficients u T and v T are complex. The phase (f>o of u T /v T corresponds to the constant shift 
of the rotation angle <j> — > (f> + (f>o in the density oscillation cos(5^>) in Eq. ©. Since the phase 4>o for the lowest-energy 
solution is isospin independent; hereafter, u T and v T are taken to be real quantities. The expectation values for this 
vacuum \^f) are 

W±k, x a ±k,x) = (b± q ,- x b±q : - x ) = u T u T , 
W±k, x b ±<i,-x) = ( b ±q,~x a ±Kx) = ±u T v T . (27) 

The normal state |0)f has v T = 1 and u T = 0, while the SSB vacuum has a finite (a^ x b± q _ x ), i.e., the finite value 
of u T v T . 

The values v T and u T are determined by minimizing the expectation value |\&), 

E = ) = H + 2E COIr , 

E co „ = J2(Ek,r - E q , T )u 2 T + 2 G^,u T v T u T ,v T ,. (28) 

Tcr x,x' 

For Z = N systems, when isospin dependences of single-particle energies Ek, T and E q ^ T are ignored, v T and u T do 
not depend on the isospin r, and the energy correction E cm from the energy Ho is 

X#X' 

= 4{(£ fc - £ 9 )u 2 + 6 5 (j,/l) uV}. (29) 
The stationary condition with respect to variations of u and v with the constraint u5u + vSv = leads to the equation 

(E k - E q )u-6g {ph) u(u 2 - v 2 ) = 0. (30) 
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For a non-zero u, u and v are solved, 

u 2 
v 2 

uv 

It turns out that, to obtain a non-zero uv with real u and v values for the SSB vacuum the following condition must 
be satisfied: 

E k -E q <-6gW. (32) 

This indicates that the SSB occurs provided that the strength —g(P h > of the attraction is large enough so as to satisfy 
the above condition. In other words, the static edge DWs can exist if the correlation energy overcomes the 

energy cost Ek — E q of a lp-lh excitation. For N = Z systems, one can regard the correlation of lp-lh as that of four 
particles, because lh corresponds to the three particles state in the particle picture. 

Let us consider the role of the proton-neutron coherence in the SSB. In the case that there is no proton-neutron 
interaction, the coupling G^fS is taken to be G^f5 = g( ph ^5 TT '(l — 5 aa >). Protons and neutrons are decoupled in the 
Hamiltonian, and the energy correction 

£corr = 4{(£ fe - E q )u 2 + 2g^U 2 V 2 } (33) 

leads to the condition for the SSB, 

E k -E q <-2gW>. (34) 

This condition is more difficult to be satisfied than Eq. (|32p . This is the reason why the proton and neutron coherent 
edge DWs can be stable, while the incoherent neutron or proton edge DW is unfavored in the oblate 7a and 3a states. 
The reason for the three times smaller interaction term, i.e., the correlation energy in Eq. (|33[) than that in Eq. (|29|) 
is that, in the particle picture, lh corresponds to the one particle state in case with no proton-neutron interaction, 
instead of the lh state corresponding to the three particles state in case with proton-neutron interactions. 

Wo also consider the further unfavored proton edge DW in a neutron-rich system discussed in the previous section. 
Protons are deeply bound in a neutron- rich system, and therefore the energy cost Ek — E q for the Ip — lh excitation 
becomes large in general. As a result, the condition Ek — E q < —2g( ph ^ becomes severe, and the proton edge DW is 
suppressed largely in neutron-rich nuclei. 

VI. SUMMARY AND OUTLOOK 

Pentagon and triangle shapes in 28 Si and 12 C were discussed in the relation with DW at the edge of the oblate 
states. In the AMD calculations, the K n = 5~ band in 28 Si and the K 71 = 3~ band in 12 C are described by the 
pentagon and triangle shapes, respectively. These negative-parity bands can be interpreted as the parity partners of 
the K n = + ground bands and they are constructed from the parity-asymmetric-intrinsic states. The pentagon and 
triangle shapes originate in the 7a and 3a cluster structures. 

We performed analysis of ideal cluster model wave functions using BB a-clustcr wave functions and also extended 
BB wave functions, and investigated the development of the pentagon and triangle shapes. It was found that the 
proton-neutron coherence is essential in development of the pentagon and triangle structures of the oblate states 
in Z = N nuclei. Without the proton-neutron coherent density oscillation, the pentagon and triangle shapes are 
suppressed. Needless to say, this is consistent with the features of light Z = N nuclei, in which cluster structures 
are favored because of a-cluster formation. In the oblate state of neutron-rich C, the triangle cluster structure is 
suppressed. 

In analysis of single-particle orbits of the AMD wave functions and BB a cluster wave functions, the pentagon and 
triangle shapes are regarded as the static one-dimensional DWs at the edge of the oblate states. The edge DWs can 
be described by nonuniform orbits with parity mixing, which give density oscillation with the wave number five and 
three at the surface of the Oku! oblate states. 




(31) 
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The static edge DWs of the oblate Z = N nuclei are understood by the spontaneous symmetric breaking (SSB) 
of rotational invariance around the symmetric axis of the oblate states. In other words, the development of the 7a 
and 3a cluster structures is interpreted as the instability of axial symmetry with respect to the pentagon and the 
triangle shapes. We introduced a simplified model and discussed the SSB for the edge DWs. The development and 
the suppression of the pentagon and triangle structures are described by the SSB inducing static edge DWs. 

In the simplified model, the OTilj oblate states are assumed to be the Hartree-Fock vacuums |0)f- The model space 
for particle and hole states are truncated so that only 4'(o > o t ±3)X Ta and 0(o,o .±2)<^t<t are active. Assuming a contact 
interaction, wc adopted the DW term Hdw as the residual interaction. For the proton-neutron coherent edge DWs in 
Z = N systems, the SSB occurs, when the condition Ek — E q < —Qg^ vh > is satisfied. If there is no coupling between 
protons and neutrons, the condition for the SSB is E k - E q < -2g (ph \ which is more severe condition than the 
proton-neutron coherent case. This means that the proton and neutron coherent edge DWs are favored, while an 
incoherent neutron or proton edge DW is unfavored. 

Considering the condition E^ — E q < —2g( ph ^ for an incoherent edge DW, we explained the reasons why the triangle 
cluster structure is suppressed in the oblate state of neutron-rich C. Since protons are deeply bound in neutron-rich 
nuclei, the level spacing of proton orbits becomes large. It increases the energy cost Ek — E q for a Ip-lh excitation, 
and hence, the correlation energy ~2g( ph ^ due to the triangle structure is not able to overcome the cost Ek — E q . 

The scenario for the suppression of the proton DW in neutron-rich systems could be extended also to infinite matter 
problems. Let us mention about the instability with respect to proton density oscillation in a neutron-rich matter, 
a symmetric nuclear matter, and a pure proton matter with the same Fermi momentum of protons ignoring the 
Coulomb force. The proton density wave should be most unfavored in the neutron-rich matter among these three 
cases, while that might be favored with the coherent neutron DW in the symmetric nuclear matter. It turns out that 
the possibility of a-cluster crystallization in the neutron-rich matter may be suspicious. Alternatively, we can say that 
the a-cluster crystallization may be suppressed in the neutron-rich matter because of the quenched effective mass of 
protons. 

In the present simplified model, in which active orbits are limited to be a small number, DW may be superior to 
BCS-type pairing in Z = N systems. We adopted the ansatz of the residual interaction H 2 = Hnw an d discuss the 
edge DWs in relation to the SSB. This ansatz may be applicable only to the case that the level density is enough low, 
active orbits are restricted in almost one dimension, and the spin-orbit force can be ignored. The oblate 12 C may 
satisfy this condition and the oblate 28 Si would do probably. However, we should comment that, in normal nuclei, 
static surface DWs may yield to the BCS pairing. The spin-orbit force may also weaken static DWs. Moreover, when 
the number of active orbits are large enough, the BCS pairing overcomes to DWs. Therefore, in heavy-mass nuclei, 
especially, in spherical nuclei, the BCS pairing can be predominant as well known. In fact, various phenomena due 
to the BCS pairing has been observed in heavy-mass nuclei and are successfully described by the BCS theory in the 
j-j coupling scheme. 



To see the relation between BB cluster wave functions and shell model wave functions, it is convenient to expand a 
BB wave function with LS'-coupling shell model wave functions which are described in terms of single-particle orbits 
in the spherical H.O. potential. For instance, an a cluster located at the origin is expressed by four nuclcons, p 
p n f, and n \. occupying the Os orbit in the H.O. potential with the frequencies ui x = us y = cu z = u) = h/mb 2 . Here 
the parameter b is related to v of cluster wave functions as v = l/2b 2 . For oblate systems, it is convenient to use the 
expression of single-particle orbits with cylinder coordinates, p = \J x 1 + y 2 , z, 0, where z is the symmetry axis. Then, 
H.O. single-particle orbits are characterized by quantum numbers n z , n p , mi- Here n z and n p are the node numbers 
with respect to z and p coordinates, m; is the eigenvalue for the z-component of the orbital angular momentum. The 
total quantum number is N = n, + 2n p + \mi\. 

The explicit forms of the H.O. single-particle orbits 0(„ zj „ pjmi ) for (n z , n p , m;) = (0, 0, ±1), (0,0, ±2), and (0,0, ±3) 
are 



Appendix A: H.O. single-particle states 




(Al) 
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Appendix B: Particle and hole representation 



In this appendix, we summarize the notations of the particle and hole representation. The creation and annihilation 
operators, c a and c a for a state \a) are defined as 

41-) = 

c a \a) = |->, 
4|a) =0, 

c a \-}=0, (Bl) 

where |— ) is the no-particle state, and a denotes the index of all degrees of freedom of the single-particle state such 
as momentum, spin, and isospin. c a and cjj satisfies {c a ,Cg} = 5 a .p, and other anticommutation relations are zero. 
To describe particle-hole excitations on the Hartree-Fock (HF) vacuum, we define the HF vacuum state as 



\0) F ^Y[d\-), (B2) 



a<F 



and the particle and hole operators as 



a a = c a for a > F, 

b a = S- a C- a for a < F, 

a a = c a for a > F, 

b a = S- a cl a for a < F. (B3) 

Here a < F and a > F means the states below and above the Fermi surface, respectively. The time reversal state of 
| a) is defined as S- a \ — a). 

For an infinite matter of spin-1/2 fermions, single-particle states can be characterized by momentum k and spin 
s z = a . In a usual convention, | — a) = \ — k, —a) for \a) = \k, a) and 

S a = (-1)^^. (B4) 

For a spherically symmetric system, single-particle states can be characterized by the quantum numbers \a) = \nlsjm,j) 
in the j-j coupling picture, and the corresponding | — a) and the phase convention are 

| - a) = \nlsj ~ m,j), 

S a = (-iy-^. (B5) 

In an axial symmetric system in the l-s coupling scheme such as the present 7a and 3a models for 28 Si and 12 C, 
we use the notation \a) = \n z n p mi<j) specified by the quantum numbers in the cylinder coordinates and adopt the 
following conventions: 

| - a) = \n z n p mi - a), 

S a = (-l)^- a - m '. (B6) 

The operator w a creates a hole carrying the z-component of angular momentum mi and the spin s z = a. 
We consider the Hamiltonian including the two-body interaction 

H = ^2{a\T\/3)c a Cf3 + - V a ,fi,y,sc a clcyC S , 

aj3 a/3y5 

V a ,p, 7 , s = \{(ap\v\it) - («0M*y>}. (B7) 

We rewrite the Hamiltonian in normal-ordered form with respect to new particle and hole operators assuming that 
the single-particle states, a, are solutions of Hartree-Fock single-particle equations, which diagonalize the Hamiltonian 
matrix 

(/3|T| 7 > = ]T [(af3\v\aS) - (a(3\v\6a)} = E p 5 0S . (B8) 

a<F 
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Then the Hamiltonian takes the form, 

H = H + H 1 +H 2 , (B9) 

with 

H = J2 <«|r|«) + i E [( a P\ v \ a P) - («0M0«>]> 

a<F a<F 0<F 

H\=y^ j E a a) a a a - ^ E a blb a , 
a>F a<F 

H ' 2 = \Y1 V<x,p,-y,8N{clc\c 5 (^), (BIO) 
a/?7<5 

where N( ) is the normal-ordered product with respect to the particle and hole operators defined before. The residual 
interaction H 2 contains the particle-particle, hole-hole, and particle-hole scattering, 

a,fi,~f,S>F 

H h h = ^ E Va^^gbib^bsby, 

a,P,~f,S<F 

H ph = 2 ^2 E ^a.-S-p^S-pS-saib^bsa^. (Bll) 

a,~/>F /3,5<F 



Appendix C: Hamiltonian of the simplified model 

We introduce a simplified model for the oblate state . In this model, <J>!^ i s assumed to be the Hartree- 

Fock (HF) vacuum |0)f, and possible particle-hole excitations are restricted within the HF single-particle states 
of |(/>(o,o ,±fe) X T a) and |<^(o,o,±g) Xra}- This means that the model space is truncated so that active orbits are only 
\4>(p,o,±k)X T cr) for particle states and \4>(o,o,±q)X T cr) for hole states with k — 3 and q = 2 (or k — 2 and q = 1 for 
the $3^"*^). We use the labels a = ±k,r<7 and a = ±q,ra for these active single-particle states and also adopt the 
notations x = Ta and — \ = t — a. In this paper, we define the particle and hole operators as 



a ±k, x ~ c ±k, x 



t 
± 

b±q,x = C T« -X- (d) 



Here, for convenience, we adopt the definition of the hole operators without the phase convention instead of Eq. (|B3|) . 

In this model, we assume a contact two-body attraction v(r) = gd(r) with g < for the residual interaction in 
the H2 term. The matrix element V a ,p,t,s is n °t zer ° only when k a + hp = k~ ( + k$ and Xa = X~t XP = Xs (or 
Xa = Xs 7^ Xp = Xy) are satisfied and calculated to be 



r<(pp) 
U xx' 


= v kx,- 


-*x'i*x.— *x' 


= g(PP){l 




g{PP) 


■f* 


-fc|<5(r)|fc,- 


-k), 




G xx' 


= v qx ,- 




= g m {1 


-<W)i 


g(hh) 


= 2^ 




-?>, 




U xx' 


= V kx ,- 


-IX 1 ,kx,-qx' 










-q\6(r)\k,- 


-3). 





(C2) 

Using the symmetry (or antisymmetry) of the matrix elements V a p~f s with respect to indexes, the Hamiltonian 
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Eq. (|B10|) in the particle-hole representation can be rewritten in the explicit form, 

H = Hq + H i + if 2 j 
H = F (0|F|0) F , 

X X 
- E E <l>X b +q,X b +<]'X ~ E E k-X b -qa b ~1'X' 



with 



H 2 = H ph + H pp + H hh , 



(C3) 



H ph = 2 J2 G 



( P h) 
x,x' 



1 +k, x b +q--x b + c i--x' a +k,x' 



+ a -k, x b -<i,-x b -<i--x' a -k,x' 
+ at +k,x b ~i,-x b -i>-x' a +k,x' 

+ a -k. X b +q,-X b +1-X' a -k,x' 
2 H G x!x' a +k a b +q-x' b +1-X' a +k,X 

x,x' 

+ al k jH q _ x ,b- g ^ x ,a- k ,x 

+ at +k,x b -q,-X ,b -^-X' a +k,X 



a t _fe, x &+ g ,_ x ^+ 9 .-x' a -fe.x 



h pp = G 



(pp) 
x,x' 



a +k,x a +k. X ' a + k -X' a +k,X 



t J 

"x,x' L" 

x,x' 

+ a -k.Y a -k.x> a -k,x' a -k, X 



a f +k ^- k)X ia-k,x' a +k,x 



(C4) 



(hh) 

x,x' 



x,x' 



b +qa b +q.X ,b +1>X' b +1>X 



b - q J-q,X 



> b -q,x' b -q,x 



& + g , x & - g , x ' & -g,x' & +<ir,X 



(C5) 



+ b +q,X b] -q,X> b +1-*' b -1-* 



(C6) 



where we omit — !> 4 and 1 — > 3 and their inverse processes. In the particle- hole interaction term H ph , the DW term 



a +fc,X°+9 -X 



fr+ g ,-x' a +fc,x' 



+ a 



t _ fc ^&L g ,- x &-g,-X' a -fc,X' 



(C7) 



may induce the edge DWs having the wave number ±(k + q), which gives the non-zero expectation value (a± fc x^±9 -x^' 
The terms of a± k x b\^ q _ x &=pq,- x 'a±fc,x' m H ph ma Y induce the exciton mode having the wave number ±1. They 
contains the spurious mode of translational motion and are of no interest in finite systems. Other terms in H ph have 
the opposite sign and they do not give coherent effects to the correlation energy. 



21 



In H pp and H , the interactions which may induce the BCS pairing are the following terms: 



ttPP 

-"bcs 



x,x' 



+ a +k, X a -k.x' a + k 'X' a -k,X 



(C8) 



Tjhh _ o \ " n {hh) 
-"bcs - 1 2^ U X,X' 
x,x' 



b +q,x b -qa' b ~9^' b +q,x 



b +qa b ~q,X' b +1^' b -q,X 



In the case of Z = N nuclei, only two types of BCS pairing, for instance, a' +k 



,t 



considered among four species of nucleons, x — P tiP 4-) n 1\ n I- This is different from the DWs induced by Huw 



and a 



+fe,nt -fc.nj. 



(C9) 
, are usually 



where four types of particle-hope combination, (a± fe x b± qi - x }, can be non zero simultaneously and they can give 
coherent effects to the correlation energy. Considering that the coupling constants, G^ ph \ G (pp \ and G (Wi) are the 
same order, the DW may be superior to the BCS-type pairing in Z = N systems in the present simplified model with 
a limited number of active orbits. We consider i?DW to be the dominant term, and adopt the ansatz of H 2 = Fdw 
and discuss the edge DWs in the Hamiltonian H = Hq + Hi + Udw in relation to the SSB. 



Appendix D: Alternative method to solve DW Hamiltonian 



In this section, we solve Eq. (|24[) in the approximation omitting the quantum fluctuation of the product of the 
particle and hole operators. It is a kind of the mean-field approaches in the field theory. We show that the same 
result as that in Sec. [V] are obtained in this approximation. 

Let us consider the H = Hq + H\ +7?dw m Eq. (|24p . By decomposing i?DW into mean fields and their fluctuations, 
we can rewrite the second and third terms H\ + i^DW as 

Hi + i^DW = H m f + i/quasi + -Hfluc, 

^^^^l^x^H-xX^r^tty) + {al k jl q ,_ x )(b- q> - x ,a- k>x ,) 



H 



quasi 



a +k, x b +q-x( b +v-x' 



a+k, X ') + ("+k,x h +9,-x) b 



+q,-x' a +k,x' 



x,x' 



+ a -k,x b ~q,~x( b -i<~x' a -k,x') + ( a -k a b ~q,-x) b ~i,-x' a -k, 



2 E« 

x,x' 



x"-i-x 

K a +k, X °+q-X 



;x' 



,t 



\ a '-k, X b -q,-X 



( a +fc,x 6 +g-x» i b +q,-x' a +k,x' ~ ( b +q,- X ' a +k,x')) 

( a -k.x b -q,-x))( b -1^X' a -k,x' - ( b -q,-X' a ~k, X ')) 



(Dl) 



where i? m f, -ffquasi, and H mi are the mean- field energy term, the sum of H\ and the interaction term between particles 
and the mean field, and the fluctuation term of the mean field, respectively. In the mean-field approximation, i?fl uc 
is assumed to be negligible, so we drop the -f/fl uc - We also assume that the ground state does not break time- reversal 
symmetry, so that the mean fields satisfy 



(b +q ,- x a+k,x) = -( a ~k,~x b -qJ- 



The Hamiltonian of the quasiparticle i? quas i can be written as 



(e Kt A 



H 



quasi 



= £(« 



l +k, x °+<?<x 



-k,x 



-q,x 



\ 



\ 






A* E q , T 
E Kt -A 
\ -A_ x E qtT ) 



( a +k,x\ 

b lq,X 

a -k, X 

\ b -qaJ 



- 2^£ , 9 , T , 



(D2) 



(D3) 



where the last term in Eq. (jD3[) comes from the anticommutation relation: b\_ q x b± qtX = —b± qyX b± q + 1. The gap is 
defined by 



Ax = 2 / ] Gx,x' (bg,-x' a k,x' 
x' 



(D4) 
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Since -ff quas ; depends on A x , the right handed side in Eq. (|D4[) also depends on A x through the expectation value; 
thus, Eq. (|D4[) can be regarded as a self-consistency equation. The Hamiltonian of the quasiparticle has the quadratic 
form, so that it can be diagonalized by the following unitary transformation or Bogoliubov transformation: 



a +k,X = l V X ~ U X a +k,X 
JhJ Vx V *x)\ b U, X ; 

a-k,x\ = ( v *x u *x\ ( a ~ fc .x 



(D5) 



with 



2Ax (D6) 



v x E Kt - E q>T + v /(S fe , T -^ r ) 2 +4|A x | 2 
which satisfies \v x \ 2 + \u x \ 2 = 1. The eigenvalues of i? qU asi corresponding to the energies of the quasiparticles are 



E Kx {A x ) = -{E Kt + E q , T + \l {E Kt - E^Y + 4|A 



2 

V 



E q . x (A x ) = -{E ktT + E qtT - yJ(E ktT - E qtT y + 4\A x \2). (D7) 
The new vacuum is defined as the state vanished by the annihilation operators, a±k, x and b± qx : 

a±fc,xl*> = b± q , x m = 0. (D8) 

The solution of Eq. (|D8|) is given by 

|*(A X )) = H(v x + u x a\ k J +kiX ) I]« - u x al k J_ Kx )\Q) F . (D9) 

x x 

This has the same form as Eq. (|25l) ; however, they arc different, because the vacuum in Eq. (|D9I) is a function of A x , 
while that in Eq. (|25j) is a function of u T (and v T ), which is a variational parameter. Their vacua coincide at the 
solutions of the self-consistent equation and the variational equation. 
The expectation value of iTq Ua si becomes 

£ quasi (A x ) = (F quasi ) = 2^(£ 9 . X (A X ) - E q>x ), (D10) 

x 

where -E quas i(A x ) < and the equality is only satisfied, when |A X | = 0. The mean-field term can be rewritten by A x 

as 

^ f =-£^^A*A x ,, (Dll) 



x,x 

where we used the explicit form of the interaction G { ^) = g& h \l - S x ^). Using Eqs. (|DT0|) and ([DTT]) . we obtain 
the correlation energy as 

2-^corr(A x ) — -^ q uasi ~l~ ^mf 

= 2E(^x(Ax) ~ E,x £ l ~^f^)- (° 12 ) 

x x' y 

In the mean-field approximation, A x is obtained by the stationary condition: 

-^co rr (A x ) = y= ~\ : = = + ^ - = 0. (D13) 



dA* x x > ^(E KT -E q , T y+±\A x \> 2g(r» 6gW 

Notice that Eq. (|D13|) is the equivalent to the consistency condition Eq. (|D4[) . which can be check by inserting 

(b q ,- x a Kx ) = u x v x . 

For Z = N systems, when E k ^ T and E q ^ T are independent of the isospin, A x is independent of x- The solution is 



(E k - E q V 

(6g) 2 



Inserting Eq. (|D14|) into Eqs. (|D6I) and (|E)12|) . one finds that the v T and u T coincide with Eq. ([3"Tj) and also the E c , 
does with Eq. (gU). 
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